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ABSTRACT 

We determine a i^-band luminosity-metallicity (K-Z) relation for dwarf irreg- 
ular galaxies, over a large range of magnitudes, —20.5 < < —13.5, using a 
combination of K photometry from either the 2-micron all sky survey (2MASS) or 
the recent study of Vaduvescu et al. (2005), and metallicities derived mainly with 
the T e method, from several different studies. We then use this newly-derived 
relation, together with published K s photometry and our new spectra of objects 
in the field of HCG 31 to discuss the nature of the possible tidal dwarf galaxies of 
this group. We catalogue a new member of HCG 31, namely "R", situated ~40 
kpc north of the group center, composed by a ring of Ha knots which coincides 
with a peak in HI. This object is a deviant point in the K-Z relation (it has 
too high metallicity for its luminosity) and its projected distance to the parent 
galaxy and large gas reservoir makes it one of the most promising tidal dwarf 
galaxy candidates of HCG 31, together with object F. The subsystems Al, E, F, 
H and R all have metallicities similar to that of the galaxies A+C and B, result 
that is expected in a scenario where those were formed from material expelled 
from the central galaxies of HCG 31. While objects Al, E and H will most 
probably fall back onto their progenitors, F and R may survive as tidal dwarf 
galaxies. We find that two galaxies of HCG 31, G and Q, have A+em spectral 
signatures, and are probably evolving toward a post-starburst phase. 

Subject headings: galaxies: dwarf — galaxies: irregular — galaxies: clusters: 
individual (HCG31) 



1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 
Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina) - Observing run ID: 
GN-2004B-Q-47. 
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1. Introduction 

It has been established for more than a decade that there is a well defined B luminosity- 
metallicity (L-Z) relation for dwarf irregular galaxies, in the sense that the higher the metal- 
licity, the higher the luminosity of the dwarf galaxy (e.g. Skillman et al. 1989; Richer & 
McCall 1995), although Hidalgo-Gamez & Olofsson (1998) claimed that the relationship was 
much weaker than previously thought. More recently the relationship was confirmed by Lee 
et al. (2003b) with a root mean square of a = 0.175 dex in log(0/H) for the same sample of 
dwarf irregular galaxies originally examined by Richer & McCall (1995), but with updated 
distance determinations and metallicity measurements. The L-Z relation is also valid for 
giant galaxies, with less scatter and with a steeper slope than for dwarf systems (Salzer et 
al. 2005). 

In the studies of dwarf systems associated with strongly interacting galaxies, this luminosity- 
metallicity relation has been used to select possible candidate tidal dwarf galaxies (TDGs) - 
newly born galaxies formed out of recycled material expelled from the parent galaxies during 
interaction. TDGs stand out from this relation, i.e., they, in general, do not follow the L-Z 
relation of dwarf irregular galaxies, but instead have an almost constant metallicity between 
1/4 and 1/3 of the solar value. Examples of B-band L-Z diagrams showing the location of 
tidal dwarf galaxies in a few interacting systems can be found in Fig. 6 of Due & Mirabel 
(1998), in Fig. 3 of Weilbacher et al. (2003) and Fig. 17 of Due et al. (2000). Nevertheless, in 
systems where star formation is widely spread and the B luminosities of the parent galaxies 
and of the TDGs may be altered, the B-band L-Z relation (used in the studies above) may 
no longer be a useful tool to select possible TDGs mainly because the B band is not a good 
tracer of the stellar mass when it is highly affected by starbursts. This is the case for the 
Hickson compact group 31, a gas- rich group, with intense star forming activity (e.g. Lopez- 
Sanchez et al. 2004), dominated by a pair of interacting dwarf galaxies, A and C, which 
clearly have their B luminosities quite affected by the light and dust involving newly born 
stars (in fact HCG 31C has Wolf-Rayet features in its spectrum; Conti 1991). In the course 
of studying the possible nature of the various sub-components of this interesting group, we 
thus felt the need of compiling from the literature a K-band L-Z relation. Such a relation 
is more useful than the B-band relation because it is less affected by a starburst, it suffers 
significantly less from absorption effects and it better characterizes the bulk of the stars (old 
component) in the group member galaxies. 

This paper is divided as follows. In Section 2, the K-band L-Z diagram for "normal" 
dwarf irregular galaxies —20.5 < Mk < —13.5 is determined from literature data. Section 3 
describes our new Gemini data, g' and r' photometry and medium-resolution spectroscopy of 
HCG 31, which are then used to determine the radial velocities, ages and metallicities of the 
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regions and to plot the K-band L-Z relation for the objects of this interacting group (the 
K photometry for the HCG 31 members come mostly from Lopez-Sanchez et al. (2004)). 
Finally, in Section 4 we discuss the fate of the TDG candidates of HCG 31. 

Distance-dependent measurements assume a distance to the group of 54.8 Mpc (derived 
by Lopez-Sanchez et al. (2004) from the Hubble law and H =75 km s _1 Mpc -1 ). We use the 
identifications for the HCG 31 objects suggested by Lopez-Sanchez et al. (2004) with one 
exception: their region E corresponds to blob E2 in this paper. We identify objects with 
single letters (e.g. objects E, F, R) and we refer to blobs which form an object with letters 
followed by numbers (e.g. blobs El and E2 compose object E). 

2. The luminosity-metallicity relation 

2.1. Previous works 

A few authors have questioned the existence of the luminosity-metallicity relation (e.g. 
Campos- Aguilar et al. 1993), specially for gas-rich objects. The B magnitude, classically 
used to investigate the L-Z relation, is known to be highly affected by strong and young 
starbursts with ages below 10 Myr, as well as by dust. Therefore, for the study of objects 
like those in HCG 31, which are known to contain young and strong starbursts, the use 
of near infrared (NIR) magnitudes should be much more robust and less affected by the 
onset of starbursts or dust absorption than the B magnitude. In the recent study of Salzer 
et al. (2005) the L-Z relation is actually explored also in the NIR regime for a sample of 
emission-line galaxies from the KISS survey. However, they present only global fits to all 
their sample galaxies, which are mostly concentrated in the high luminosity part of the 
diagrams (M^ < —21), whereas the small numbers of dwarf galaxies in their sample appear 
to follow a relation with a shallower slope. Recent efforts directed to test the existence of an 
L-Z relation for dwarf irregular galaxies in the H-band, by following an approach conceived 
to minimize the effects of uncertainties in distance determinations (Saviane et al. 2005) 
have led to encouraging results, although still based on a small number of galaxies. Here 
we investigate the L-Z relation in the K band by compiling from the literature a sample 
of nearby dwarf irregular galaxies selected for spanning a wide range in luminosity, having 
oxygen abundances obtained following an homogeneous method (as far as possible), and 
having reliable distance determinations. 
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2.2. Data used in our compilation 

In Table 1 we list absolute K s magnitudes and oxygen abundances for a number of 
nearby irregular galaxies. The sample of 29 galaxies was taken either from the recent work 
of Vaduvescu et al. (2005), where new NIR measurements of dwarf irregular galaxies are 
presented, or from the classic work of Richer & McCall (1995) and includes two additional 
galaxies NGC 1705 and NGC 1156 (taken from Hunter & Elmegreen 2004). It contains only 
nearby (distance < 8 Mpc) irregular galaxies, for which distances can either be determined 
from the brightest-stars (bs) method, tip-of-the-RGB techniques (rgb) or Cepheids (cep). 
For one of the galaxies, the available distance was determined from the Tully-Fisher relation 
(tf). All distances were obtained from the compilation of Karachentsev et al. (2004), except 
for the distance to IC 4662, which was taken from Lee et al. (2003b). 3 NIR magnitudes 
for the galaxies not included in the sample of Vaduvescu et al. (2005) are from the 2MASS 
extended source catalog (Jarrett et al. 2000), except for NGC 5408 whose K s magnitude is 
taken from Noeske et al. (2003). 

The values of foreground Galactic extinction in the i^-band (A^) for the sample galaxies 
are listed in column 3 of Table 1. These are taken from NED and were determined following 
Schlegel et al. (1998). Exceptions are galaxies NGC 1569 and IC 10 that are located at 
low Galactic latitudes, where extinction values derived from the dust maps of Schlegel et 
al. (1998) are not accurate. For these galaxies we adopted the Galactic extinction values 
determined by Origlia et al. (2001) and Richer et al. (2001), respectively. 

The oxygen abundances, compiled from the literature, were determined in two ways. 
Column 6 of Table 1 shows determinations based on electronic temperatures measured from 
the [O III] A 4363 line (T e -method), and column 7 shows determinations from the empirical 
calibration of Pilyugin (2001a,b, P-method). An exception is NGC 3738 for which only 
a metallicity value determined with the method of Edmunds & Pagel (1984) was found. 
In a few cases, P-method abundances were not available in the literature, but we could 
use published values of emission-line intensity ratios (corrected for extinction) to calculate 
oxygen abundances following Pilyugin (2001a,b). In particular we used line intensities from 
Lee et al. (2003b) for IC 10 and NGC 1560, from Thuan & Izotov (2005) for Mrk 209, from 
Kennicutt & Skillman (2001) for NGC 4789A, from Masegosa et al. (1991) and Miller & 



3 We did not include the following galaxies from the list of Vaduvescu et al. (2005), because we were not 
able to locate the corresponding measurements for their mctallicities in the literature: Cas 1, Mbl, Orion 
Dwarf, UGC 4115, UGC 4998, UGC 5692, UGC 5848, UGC 8508, UGC 5979, NGC 3741, NGC 4163, NGC 
4190, and Holmberg IV. On the other hand, the dwarf irregular galaxies from the list of Richer & McCall 
(1995) Sextans A, Sextans B, LMC, SMC, WLM, Leo A, IC 1613, and NGC 2366 were not included because 
no NIR total magnitudes are available for them. 
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Hodge (1996) for IC 2574 and DDO 50, and from Masegosa et al. (1994) for NGC 5408. 

2.3. The K-band L-Z diagram for dwarf irregulars 

Abundances against extinction corrected K s magnitudes are plotted in Fig. 1 (Table 1, 
values of column 6 plotted against those of column 2, after correction with extinction values 
of column 3). Data for all sample galaxies are used with T e -method abundances, when 
available (i.e. for all but two galaxies, see Table 1, column 6). A non-weighted least-squares 
fit to the data gives the following result: 

12 + log(0/H) = (-0.14 ± 0.02) x M Ks + (5.55 ± 0.26), (1) 

with a a of 0.15 dex in log(0/H). 

To establish whether there is a correlation in the data, we calculated the Spearman's 
rank correlation coefficient, a non-parametric measure of correlation, and found a value of 
—0.88. This indicates the presence of a good anti-correlation, whose level of significance can 
be inferred by comparison with published tables (Wall & Jenkins 2003). Taking into account 
the number of L-Z pairs in the sample, the above correlation coefficient indicates that the 
hypothesis that the variables are unrelated is rejected at the 0.1 per cent level of significance. 
Residuals of the fit are shown in the lower panel of Fig. 1. 

A good correlation is also present if we plotted the subsample of 25 dwarf galaxies that 
have abundances determined through the P-method (not plotted here). Also in this case the 
correlation is good 4 , with a Spearman-rank correlation coefficient of —0.81, and same level 
of significance as above. The result of the non-weighted least-squares fit to the data is not 
significantly different from the one above: 

12 + log(0/H) = (-0.13 ± 0.02) x M Ks + (5.73 ± 0.34), (2) 

with a a of 0.16 dex. 



3. Observations and results for HCG 31 

We have used relation (1) above, combined with available K-band photometry, from 
(Lopez-Sanchez et al. 2004), and our new medium- resolution spectroscopy of the objects in 



4 We excluded from the fit galaxy NGC 6822, whose metallicity was determined with the P-method, and 
which deviates considerably from the L-Z relation (Pilyugin 2001b). 
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HCG 31, to study the nature of the members of this group, as deduced from their location 
in the K-band L-Z relation for normal dwarf galaxies. Our new data are described in the 
following. 

3.1. New data - Observations 

New imaging and multi-slit spectroscopic observations of HCG 31 were done with the 
GMOS instrument, mounted on the Gemini North telescope, on August 29 and September 
21 of 2003, respectively. 

The imaging consisted of 5 x 180 s exposures in the r' band, and 5 x 240 s exposures 
in the g' band. The filters are from the SDSS system (Fukugita et al. 1996). The typical 
FWHM for point sources was ~ C/75 in all images. The observations were performed in 
photometric conditions. Fig. 2 displays the r' image of HCG 31. Fig. 3 displays zoomed 
pannels of selected regions, with 20 surface brightness contour levels logarithmically spaced 
from 17.8 to 24.0 mag arcsec -2 . The positions of the spectroscopic slits are also indicated. 

Standard reduction steps were performed with the Gemini package GMOS. After flat- 
fielding and cleaning from cosmic-ray events, the final frames were analyzed with the program 
SExtractor (Bertin and Arnouts 1996). 

The calibration to the standard SDSS system was made with the general zero points and 
extinction coefficients provided by the Gemini observatory 5 . The accuracy of the calibration 
is claimed to be within 5% to 8%. 

Three multi-slit exposures of 960 seconds each were obtained through a mask with 1.0" 
slits, using the R400 grating, for a final resolution of 6.0-6.5 A, covering approximately 
the range 4000 - 8000 A. Three additional multi-slit exposures of 1200 seconds each were 
obtained through a mask with 1.0" slits, using the B600 grating, for a final resolution of 4.5 
A, covering approximately the range 3750 - 6600 A. The typical FWHM for point sources, 
measured on images taken for identification of the field, was ~ 0"6. Only a few of the 
observations were performed in photometric conditions. 

Standard procedures were used to reduce the multi-slit spectra using tasks within the 
Gemini IRAF 6 package. Wavelength calibration was done using Cu-Ar comparison-lamp 



5 www . gemini . edu/ sciops/ instruments/gmos/gmosPhotStandards . html 

6 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
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exposures before and after the exposure on the target. Flux calibration was done using 
spectroscopic standard stars obtained in the same night of the observations. The blue and 
red spectra were glued together, after flux calibration, and are shown in Fig. 4. 

3.2. Measured properties of the HCG 31 objects 

In columns 2 and 3 of Table 2 we list the coordinates of all objects in the HCG 31 group 
identified in the r' image. We include a new object which we named "R", at RA = 05 h 01 m 
34 s and DEC = -04°12'57" (JD2000), located 2.5 arcmin north of the group central object 
C, about 1 arcmin northwest of object Q (see Figs. 2 and 3). As discussed in Section 4, this 
may be one of the best candidates for a tidal dwarf galaxy in HCG 31. We only obtained a 
spectrum for one of the blobs which constitute region R, namely Rl. 

The remaining columns of Table 2 list measurements made by us and other authors on 
the properties of the objects in the group, such as magnitudes, Ha-luminosities, velocities, 
metallicities and colours of the group members. Column 4 lists the K absolute magnitudes, 
computed from the K apparent magnitudes measured in Lopez-Sanchez et al. (2004), for an 
adopted distance to the group of 54.8 Mpc (values for regions H, Q, and R were measured by 
us, see details below), column 5 lists the logarithm of the luminosity in Ha, when available, 
from Lopez-Sanchez et al. (2004) or from this work (in the latter case they are not corrected 
for light loss from the slit, and are therefore lower limits), columns 6 and 7 list respectively 
the heliocentric radial velocities of the objects (with errors) from this work and from the 
H I velocities obtained by Verdes-Montenegro et al. (2005), when available. The next four 
columns of Table 2 list four different determinations of 12 + log (O/H): using the T e method 
and the N2 estimator (results from Lopez-Sanchez et al. 2004), in columns 8 and 9, and 
using the N2 and 03N2 estimators, derived from our own data, in columns 10 and 11 (see 
the definition of these metallicity estimators in the following subsection). Columns 12 and 
13 list the measured line ratios used in the determination of the metallicities of columns 10 
and 11. Equivalent widths of Ha are listed in column 14 and g'-r' colours (measured within 
an aperture of 2") are in the last column of Table 2. Details are presented below. 

Radial velocities: 

The spectra of all objects marked in Fig. 3 (with exception of A, Al and A2) are shown 
in Fig. 4. All spectra have emission lines. The heliocentric velocities of the observed objects, 
derived from the redshifts of the brightest lines, are listed in column 6 of Table 2 (the errors 



Science Foundation. 
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are the rms of the individual line measurements). The H I velocities at the location of each 
corresponding object (from Verdes-Montenegro et al. 2005) are given in column 7. We note 
that the velocities of all optical regions studied here coincide with the H I velocities from the 
channel maps within the errors (except for region C), suggesting a physical association of the 
objects with the H I clouds. The velocities are in agreement with those derived from long-slit 
spectroscopy by Rubin et al. (1990) and Hickson et al. (1992), for the objects measured in 
common. 

Metallicities: 

Besides the metallicity estimates obtained by Lopez-Sanchez et al. (2004), not available 
for all objects (see columns 8 and 9 of Table 2), we have computed values for 12 + log (O/H) 
using two metallicity indicators, from our own measurements of the line ratios. The results 
are listed in columns 10 and 11 of Table 2. 

Our first estimate was obtained with the N2 calibrator, following Pettini & Pagel (2004), 
which is defined as the logarithm of the [NII]A6584/Ho; ratio. The resulting metallicities and 
the measured ratios are listed in columns 10 and 13 of Table 2 respectively. Our values in 
column 10 can be directly compared with those of column 9, derived by Lopez-Sanchez et 
al. (2004). As can be noted, these completely independent measurements are very similar, 
for the objects measured in common. 

Our second estimate was made using the 03N2 index (Pettini & Pagel 2004), based on 
the logarithm of the ([OIII]A5007/H/3)/([NII]A6584/Ha) ratio. The resulting metallicities 
and measured ratios are listed in columns 11 to 13 of Table 2. 

The rms scatter in the calibration of these estimators is ~ 0.25 - 0.4, which is larger 
than the internal errors; we then assume that 0.3 is the one-sigma error of our measurements. 
The metallicities inferred for the objects are very similar within the errors, the mean value 
being 12+log(0/H) = 8.3 for our estimates (with either method). 

We point out that the metallicity obtained for object Q is very uncertain given that 
the weak emission lines are superposed onto strong absorption lines, hampering a reliable 
determination of the metallicity for this object. We include this object in the tables but do 
not plot its magnitude/metallicity in the K-Z relation of Fig 1. 

Equivalent widths and ages of objects E, F, H and R: 

For objects El, Fl, F2, F3, HI, H2 and Rl, where hardly any continuum is detected (see 
Fig. 4), we assume that the objects are excited by young stars formed in an instantaneous 
starburst and use Starburst99 (Leitherer et al. 1999), with solar metallicity and Salpeter 
IMF, to estimate ages (obtained from the equivalent widths of Ha, as given in column 14 of 
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Table 2). 

From a comparison between the observed Ha equivalent widths and those produced by 
Starburst99, we find ages around 3 Myr for all regions but F3 (the latter has an approximate 
age of 6 Myr, while the youngest one, Rl, has an age of 2.6 Myr). 

Spectroscopic properties of Q and G: 

The two galaxies Q and G, besides exhibiting H II-region-like emission lines, also have 
strong Balmer absorption lines, typical of spectra dominated by A- and early F-type stars, 
with B.5 equivalent widths significantly larger than that of normal spiral galaxies. The 
equivalent widths (EW) of the [OII]A3727 emission line and the B.5 absorption line are often 
used to classify galaxy spectra on the basis of their current/past star formation episodes 
(e.g. Poggianti et al. 1999; Poggianti & Wu 2000; Mateus & Sodre 2004). Galaxy G, with its 
EW([0 II]) = ~ 20 A and EW(FPj) ~ 8A, falls in the category of the so-called e(a) galaxies 
(Poggianti et al. 1999), or A+em galaxies, according to the notation of Balogh et al. (1999), 
with still considerable ongoing star formation. The spectrum of galaxy Q does not include 
the [OII]A3727 line, however its weak Ha emission line indicates that a low level of current 
star formation is still present. Its spectrum is dominated by A-type stars, as indicated by the 
particularly strong Balmer absorption lines, EW(H5) ~ 13 A, a rather extreme and unusual 
value. Also this galaxy can then be considered of e(a) type. 

Aperture (g'~ r') for all objects and K magnitudes of H, Q and R: 

Aperture magnitudes in the g' and r' bands were obtained for all studied objects, within 
an aperture of 2 arcsec, using the task phot in IRAF. 

K magnitudes for objects H, Q and R, which were not given in Lopez-Sanchez et al. 
(2004) were derived by us using the program Sextractor (Bertin and Arnouts (1996), param- 
eter magbest). We measured these three objects in unpublished J images of HCG 31 and 
assumed a J-K colour of 0.345 for all the objects. We used J-band images from the archive 
of the New Technology Telescope of the European Southern Observatory, obtained with the 
instrument SOFI. Nine images of HCG 31 obtained on Nov 3rd/2001 were retrieved, of which 
only five contained objects Q and R (object H was visible in all nine). Images in the K band 
were also available in the archive but they did not go deep enough to allow measurements of 
the objects. Calibration of the instrumental photometry was done by using stars common 
in our fields and in the 2MASS images of the group. Our final magbest J magnitudes were 
18.6±0.1, 14.80±0.04 and 19.7±0.5 for objects H, Q and R respectively. These values were 
corrected for the foreground Galactic extinction given in NED, Aj = 0.046 mag and trans- 
formed into K magnitudes using a J-K colour of 0.345 (which is an average of the J-K colour 
of the other members of the group). The corresponding K absolute magnitudes for the three 



- 11 - 



objects, assuming a distance to the group of 54.8 Mpc, are listed in column 4 of Table 2. 

3.3. The K-Z relation for the TDG candidates of HCG 31 

We overplot in Fig. 1 the values for the metallicity and K magnitude (values in columns 
8, 9 or 10 against those of column 4 of Table 2) for the components of HCG 31. The data 
on the metallicities come either from (Lopez- Sanchez et al. 2004) or from this work. Those 
from this work (for objects H and R only) were derived through the N2 calibrator (column 
10 of Table 2). The only value plotted from column 9 was that for the region Al, for which 
no other measurement of metallicity is available. 

In Fig. 14 of Lopez-Sanchez et al. (2004), which plots the B-metallicity relation for the 
HCG 31 members, galaxies B, C and G are more than 2 magnitudes off the line for the B-Z 
relation for normal dwarf galaxies. In our corresponding Fig. 1, using K magnitudes instead, 
these same objects are closer to the best line that fits the relation for normal dwarfs, while 
still in the high-luminosity side. On the other hand, for the fainter objects, only objects H 
and R stand off the correlation, although also Al, E2 and F1+F2 show a larger metallicity 
and/or a fainter magnitude than the best-line fit. We suggest that A, B, C, G and Q are 
galaxy members of HCG 31 while the fainter objects are either tidal debris or tidal dwarf 
galaxies of the group formed due to the interaction. This is further discussed in section 4. 

4. Discussion 

4.1. General Properties of HCG 31 members 

Table 3 summarizes all of the properties of the HCG 31 members either from this work 
or gathered from the literature. This table may be important for future modelling of the 
group, for comparison between simulations and observations. As can be noted, HCG 31 
has been observed almost in all wavelengths. The optical morphology of the galaxies in the 
group is very disturbed, as is also the optical and HI kinematics. 

A point which is worth noting is the large difference between the equivalent widths of 
Ra of object F3 (EW = 74 A ), as compared to those for Fl and F2 (1508 A and 1010 
A respectively). This, in turn, leads to a large difference in the derived ages of F1+F2 (3 
Myrs) and F3 (6 Myrs). The colours g'-r' of the three blobs are very similar and are the 
bluest colours observed for any object in this group (see last column of Table 2). Amram 
et al. 2004 showed that the rotation curve is flat through Fl and F2 and the velocity of F3 
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is ~ 40 km s -1 higher than that for Fl and F2. The values for the metallicities are similar 
for the three blobs (within the errors) although F3 tends to have higher metallicity than the 
other two. F3 is also the lowest surface brightness object of the three. One might think 
that the different equivalent widths and the (small) discontinuity in the kinematics could be 
hints that F1+F2 is a distinct object from F3, but we not not believe this is the case, given 
the detailed analysis of the kinematics of the region: we have revisited the Fabry-Perot data 
cube of HCG 31 (see Amram et al. 2004 for a description of the data) and have seen that 
there clearly exists a continuity in the velocity field between F1+F2 and F3. The situation 
for object E is also similar: although the equivalent width of Ha for El (740 A ) is much 
larger than that for E2 (34 A ), analysis of the kinematics done in Amram et al. (2004) 
shows that E1+E2 form one single object. 

One new idea put forward in section 3.2 and also listed in Table 3 may deserve some 
discussion. Although the optical spectra taken in this study show mainly HH-region-like 
spectra for the HCG 31 regions (see Fig. 4), for two of the galaxies, G and Q, the spectra 
are typical of e(a)-type galaxies. The properties of the e(a) class of galaxies have been 
interpreted by Poggianti et al. (1999) as a possible indication of dusty starbursts. However, 
in the spectra of Q and G there are no indications of particularly high values of internal 
extinction affecting the emission lines, since the observed Ha/H/3 ratios are moderate (~ 
3.4 for Q and 4.0 for G). Hence, the interpratation of these two spectra remains uncertain. 
Nevertheless, we suggest that galaxies Q and G have started reducing their star formation 
rate and are evolving towards a post-starburst phase. 



4.2. The tidal tails and the TDG candidates of HCG 31 

HCG 31 is completely embedded within an H I envelop and it has two main tails: 1) the 
southern tail, which contains objects E, H, F and G, is a narrow and linear optical and Ha 
tidal tail, starting from galaxy C towards the southeast, ending with galaxy F or perhaps 
going even further (to the region of diffuse emission situated between F and G). This tail has 
been the subject of several previous studies (e.g. Amram et al. 2004) and 2) the northern 
tail, including objects Al, Q and R (see Fig. 2). The base of this tail, close to galaxy A, 
shows an open configuration, suggesting that the material has moved from its original plane. 
Except for the H I study of Verdes-Montenegro et al. (2005), this northern tail has never 
been studied before. 

N-body simulations for compact groups show that stellar tidal tails are transient features 
that can be easily destroyed due to multiple interactions (e.g. Athanassoula et al. 1997; 
Barnes 1985). Gaseous tails may also have similar fates. The frequency of occurrence and 
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length of the tidal tails in galaxy mergers are strong functions of the encounter geometry 
and the merger phase. Galaxies in the pre-merger phase, where the two galaxies are still 
distinct but have gone through the first encounter, are expected to have well developed tidal 
features, as seen for HCG 31. As the final merger takes place, the tidal tails gradually 
disappear, with the material in the tails being accreted back onto the remnant or escaping 
the system altogether, eventually forming TDGs (e.g. Mihos et al. 1998; Hibbard & van 
Gorkom 1996). 

In Fig. 1 we plotted the location of the components of HCG 31 in the K-Z diagram in 
an attempt to identify good TDG candidates in this group. The objects display a range of 
luminosities of at least 6 magnitudes and a range of oxygen abundances of ~ 0.5 dex. 

While in the B-Z relation plotted in Fig. 14 of Lopez-Sanchez et al. (2004) objects B, 
C and G are more than two magnitudes off the relation, in the corresponding K-Z relation 
these objects are closer to the best line for normal dwarfs, although they are still in the 
high-luminosity side of the relation. Although objects Al, E2 and F1/F2 are located in the 
low-luminosity side of the relation, they do not stand out, perhaps because even their K 
magnitudes could be brightened by the strong star bursts present in these objects (it is well 
known that at least F has a low old-stellar population content, Johnson & Conti 2000). On 
the other hand, the low-luminosity objects H and R are completely off the K-Z relation. We 
note that for galaxy Q the metallicity determination has a very large error given the weak 
emission lines in its spectrum - this is then not plotted in Fig. 1. 

We suggest, based on the velocities, positions in the L-Z relation, on their morphologies 
(our Figs. 2 and 3) and on their internal kinematics (Amram et al. 2004), that A, B, C, G 
and Q are galaxy members of HCG 31. On the other hand, based on the arguments below, 
we suggest that the lower-luminosity objects are either tidal debris or tidal dwarf galaxies 
of the group, formed as a consequence of the interaction. The more difficult task is then to 
decide which of the lower luminosity objects Al, E, F, H, R are tidal dwarf candidates or 
merely tidal debris. 

Iglesias-Paramo & Vilchez (2001) devised a scheme to pick out tidal dwarf galaxy can- 
didates in compact groups, based on their projected distances to the nuclei of the parent 
galaxies (at least 2 R25) and their Ha luminosities (which should be greater than 10 38 ergs 
s _1 ). We note that HCG 31 was present in the sample of Iglesias-Paramo & Vilchez (2001) 
and our regions E2, HI, H2, Fl, F2, and F3 correspond to their regions c, d+e, f, g, h, and 
i. The last three, composing object F, were among their final list of good candidates for 
TDGs. Following the same criteria, we classify one additional region as a tidal dwarf galaxy 
candidate, namely, region R (not known at the time of that study). 
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Verdes-Montenegro et al. (2005) also noted that object F is a good candidate for a tidal 
dwarf in formation. HCG 31 Fl coincides with a peak H I column density of 3 x 10 21 atoms 
cm -2 . In addition, the three blobs, Fl, F2 and F3 are the bluest objects of the group (see 
last column of Table 2). No underlying old stellar population has been detected for this 
object (Johnson & Conti 2000) and it has a metallicity similar to that of the central galaxies 
of HCG 31 (A+C and B). Although Amram et al. (2004) measured no rotation for F1+F2, 
which could suggest that it may not turn into an independent object, recent simulations of 
tidal dwarf galaxy formation (Wetzstein et al. 2005) conclude that these objects are expected 
to be non-rotators or very slow rotators. 

Object R is a second excellent candidate tidal dwarf galaxy in the HCG 31 group. It 
coincides with an H I cloud of column density of 10 21 atoms cm -2 , which is one of the highest 
seen among all the tidal filaments of HCG 31. Optically it is formed by an approximately 
round distribution of faint Ha knots. We obtained a spectrum for one of these knots (Rl), 
and we confirm that it is at the redshift of the group. Different from the intergalactic HII 
regions observed by Mendes de Oliveira et al. (2004), these knots are located in the peak of 
the HI distribution and seem to be a chain of linked star forming regions. 

Objects F and R are then good examples of tidal fragments, well separated from their 
progenitors, coincident with peaks of the HI distribution, that could evolve to become tidal 
dwarf galaxies. Objects Al, E2 and H are likely also tidal fragments, but much closer to 
their progenitors, and may be falling back to the main central object of the group. In fact, 
E has two components in counter- rotation (Amram et al. 2004), suggesting that some of the 
material may be already returning to the parent galaxy. 

Our hypothesis that regions F and R are good TDG candidates is also supported by the 
comparison of their properties with those of tidal dwarf galaxies that have been identified in 
other systems, for example those studied by Iglesias-Paramo & Vflchez (2001), Weilbacher 
et al. (2003), and by Temporin et al. (2003). The TDGs identified by Iglesias-Paramo & 
Vilchez (2001) in a subsample of HCGs have Ha luminosities comparable to those of the 
TDG candidates studied in this work. Similar results concerning Ha luminosities of TDG 
candidates were obtained by Weilbacher et al. (2003), who found an average Ha luminosity of 
2.2 x 10 39 ergs s" 1 for TDGs identified in a sample of ten interacting systems, with the most 
luminous knots having values of order 10 40 ergs s -1 . Other candidate TDGs were identified 
by Temporin et al. (2003) in the tidal tail of the compact group CG J1720-67.8, whose 
properties seem to indicate an evolutionary stage similar to that of HCG 31 (Temporin & 
Fritze-v. Alvensleben 2005). The estimated masses of the TDGs in this system are of order of 
2 x 10 7 M and burst ages for most objects range from ~ 6 to 8 Myr. Ha luminosities from 
integral field spectroscopy are of order of 10 40 erg s _1 , with the brightest TDG candidates 
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having L(Ha) = 4 x 10 40 erg s -1 , without correction for internal extinction (Temporin 
et al. 2005). The above ages are comparable to that obtained for object F3 (c.f., Section 
3.2) and the Ha luminosities are of the same order of magnitude of those of our brightest 
components, listed in column 5 of Table 2. Additionally, we note that TDGs are usually 
associated with H I density peaks (e.g. TDGs in M81 and in NGC 5291, Makarova et al. 
2002; Due & Mirabel 1998, respectively), as it is the case for objects F and R (see fig. 3 of 
Verdes-Montenegro et al. 2005). 

Regarding the metallicities, when we place the HCG 31 objects in the diagram of Fig. 1, 
we note that for the lowest luminosity objects (H and R), the metallicites we measure are 
much higher than those for normal dwarf irregulars of similar luminosities, although F is 
not far from the best line for dwarf galaxies. We suggest that mostly the gas and some 
stars which today form regions F and R were earlier in galaxies A+C or B and were torned 
out by the interaction. Then in more recent times (a few million years ago) the young star 
complexes in F and R were formed through compression of the intergalactic H I gas by galaxy 
collisions. In fact, the young star complexes in R and F would have had to be formed in 
situ and not be ejected from the central galaxy because the time needed for those to move 
even from the outskirts of the central interacting pair to their current location, at a typical 
ejection velocity (due to dynamical interactions) of 200-300 km s" 1 would be much longer 
than the age of their massive stars (as determined in Section 3.2). 

N-body simulations of galaxy collisions have shown that tidal dwarf galaxies often form 
in the tails of major merger galaxies. The timing of each star-forming burst along the tails 
is strongly determined by the orbital orientation (prograde or retrograde) and the internal 
structure of the merging galaxies (bulge or bulgeless). Recently, through high resolution 
N-body/SPH simulations, Wetzstein et al. (2005) have shown that bound stellar objects can 
only form in the tidal arms of interacting disk galaxies, if they have a sufficiently massive 
and/or extended gas component. To our knowledge, no simulation including more than three 
galaxies, as is the case for compact groups, has ever been done, to attempt to reproduce an 
specific observational configuration. Simulations of HCG 31 are highly needed to understand 
the evolution of this complex interacting system and the formation of new objects. 
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Fig. 1. — K s L-Z relation for a sample of 29 dwarf irregular galaxies (filled triangles) out of 
which 27 have oxygen abundances determined with the T e -method and two have abundances 
determined from bright emission lines (Table 1). Overplotted are the components of HCG 
31, indicated by the relevant labels. Residuals of the fit are shown in the lower panel. 
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Fig. 2. — Image in the r' band of HCG 31 obtained with Gemini+GMOS, in a logarithmic 
gray-scale representation. North is up and East is to the left. The horizontal bar indicates 
1 arcmin. The individual group members are identified. All regions marked on the figure 
belong to the group (i.e. they are at the same mean redshift of 0.013). 
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Fig. 3. — Zoomed r' images of various regions of Fig. 2, to show details of the morphology, 
identification of the components and placement of the spectroscopic slits. The horizontal 
bar indicates a distance of 10 arcsec. The 20 surface brightness contour levels are spaced 
logarithmically and go from 17.8 to 24.0 mag arcsec" 2 . The location and orientation of the 
spectroscopic slits are also shown. 
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Table 1. K„ Magnitudes and Metallicities of Dwarf Irregular Galaxies 
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3,4 


Mrk 209 a 


-16.27 


0.005 


28.78 


bs 


7.82 


7.76 
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14 
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10, this work 
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rgb 
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15, this work 
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-15.16 
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UGC 4483 
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rgb 
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7.47 


5,4 


UGC 5423 


-16.53 


0.029 
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Note. - Column 1: Identification of the object. Column 2: K s magnitudes either from 
Vaduvescu et al. (2005) or from 2MASS, before correction for Galactic extinction. Column 
3: Galactic extinction from NED or from Origlia et al. (2001, NGC 1569) and Richer et al. 
(2001, IC 10). Column 4: Distance moduli (D. M.) from Karachentsev et al. (2004). Column 5: 
Methods used for distance determination: 'brightest stars' (bs), 'tip of RGB' (rgb), and Tully- 
Fisher relation (tf). The D. M. for IC 4662 is from Lee et al. (2003b). Columns 6 and 7: Oxygen 
abundances derived with the T e -method and the P-method, respectively. When abundances were 
present in the literature for more than one H II region of a given galaxy, the average value was 
considered. Column 8: References for the metal abundances; when two references are indicated, 
the first is for the T e -method and the second for the P-method determination. 

a These galaxies have been classified as blue compact dwarfs by some authors: Richer et al. 
(2001, IC 10); Hunter & Hoffman (1999, Mrk 209); Noeske et al. (2003, Mrk 178, IC 4462, 
NGC 5408); Gil de Paz et al. (2003, NGC 1705) 

b Lee et al. (2003b) derived for NGC 1560 a lower limit 12+log(0/H) > 7.97 

c Metallicity derived with the empirical method of Edmunds & Pagel (1984) 

d As noted by Lee et al. (2003a) the P-method value for this galaxy is uncertain, although in 
agreement with other methods of abundance determination. They adopt as best value an average 
of the values given by different methods, 12+ log(0/H) = 8.61, and note that this galaxy has an 
anomalously high metallicity for its luminosity. 

e K s magnitude from Noeske et al. (2003) 

f Pilyugin (2001b) note that the H II regions of IC 4662 seem to belong to the transition zone 
of the R23 - O/H diagram, hence, in strict sense, the P-method equation for low metallicity H II 
regions should not be used. A similar argument holds for NGC 5408. 

References. — (1) Pilyugin (2001a); (2) Mateo (1998); (3) Hidalgo-Gamez & Olofsson (1998); 
(4) Pilyugin (2001b); (5) Skillman et al. (1994); (6) Miller & Hodge (1996); (7) Lee et al. (2003b); 
(8) Hunter & Hoffman (1999); (9) Thuan & Izotov (2005); (10) Kennicutt & Skillman (2001); 
(11) Hidalgo-Gamez & Olofsson (2002); (12) Lee et al. (2003a); (13) Hidalgo-Gamez et al. (2003) 
; (14) Skillman et al. (1989); (15) Masegosa et al. (1994); (16) Lee & Skillman (2004); (17) Richer 
& McCall (1995); (18) Vigroux et al. (1987); (19) Masegosa et al. (1991) 
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0.42 


El 


5:01:37.23 


-4:15:48.9 




40.05-39.15 


3967+ 10 




8.13 


8.36 


8.31 


8.30 


2.095 


0.092 


740 


0.59 


E2 


5:01:37.48 


-4:15:57.5 


-16.28 


38.90 c 


4009± 14 








8.30 


8.25 


2.770 


0.088 


34 


0.42 


Fl 


5:01:39.71 


-4:16:22.2 


-15.97 


40.12-40.01 


3969± 9 


3968 


8.07 


8.06 


8.00 


7.98 


5.850 


0.026 


1508 


0.16 


F2 


5:01:40.15 


-4:16:27.5 


-14.67 


39.67-39.74 


3956± 9 




8.03 


8.02 


7.98 


7.99 


5.027 


0.024 


1010 


0.19 


F3 


5:01:40.62 


-4:16:36.5 




37.87 c 


3990± 90 








8.13 


8.29 


1.033 


0.044 


74 


0.13 


G 


5:01:44.01 


-4:17:19.5 


-19.67 


40.60-39.16 


3997± 9 


4005 


8.15 


8.41 


8.44 


8.35 


2.325 


0.155 


22 


0.53 


HI 


5:01:38.20 


-4:16:07.2 


-15.5° 


38.17-38.64 


4050± 29 






8.41 


8.33 


8.52 


0.459 


0.102 


600 


0.46 


H2 


5:01:38.72 


-4:16:13.5 




38.59 c 


4066± 23 








8.24 


8.23 


2.444 


0.069 


556 


0.42 


Q 


5:01:38.30 


-4:13:20.9 


-19.3 


38.11 c 


4037+125 


4090 






8.51: 


8.56: 


0.694 


0.207 


6 


0.56 


Rl 


5:01:34.33 


-4:12:56.7 


-14.3* 


38.12 c 


4139±300 








8.19 


8.30 


1.231 


0.057 


1672 


0.44 



Note. — Columns 6 to 9 list four independent measurements of 12+ log(0/H), column 6 from Lopez-Sanchez ct al. (2004), using the T e method (row 3 
of their Table 5), column 7 using the NII/Ha ratio from Denicolo et al. (2002) (row 5 of Table 5 of Lopez-Sanchez ct al. (2004)). Column 8 lists results of 
our measurements using the NII/Hct ratio (listed in column 11) and column 9 lists those using the N203 ratio (column 10 over column 11), both indicators 
denned in Pettini & Pagel (2004). 

a From Table 2 of Lopez-Sanchez ct al. (2004), for an assumed distance of 54.8 Mpc, except for the measurements for H, Q, and R, that were estimated by 
us from an unpublished J-band image, corrected for the foreground Galactic extinction given in NED, Aj=0.046 mag, and adopting J-K=0.345. 

b From Lopez-Sanchez et al. (2004), spectroscopic measurement with corrections 

c Our measurements; these arc lower limits because the luminosities were not corrected for light loss from the slit. 
d From Table 5 of Lopez-Sanchez et al. (2004) 

e This value for the magnitude corresponds to the whole region H 
f This value for the magnitude corresponds to the whole region R 



Table 3. History and Properties of HCG 31 members from Literature Data 



Nature Properties References 



HCG 31 AC 



Galaxies in pre- 


Underlying old stellar population 


1 


merging phase 


Wolf-Rayet galaxy 


3 




SB nucleus moving to post-SB phase 


4 




Major interaction occurred 400 Myr ago 


1 




Age of central burst 4-5 Myr 


1 




Young massive clusters < 10 Myr 


1, 2 




Bulk of 7.7 and 14.3 /xm emission 


4 




Low 14.3 /im to 7.7 flux ratio 


4 




Contain bulk of dust of the system 


4 




Peak of CO and FIR in overlapping region of A and C 


7 




Double kinematic components 


5 




Double nucleus 


6, 2 


HCG 31 B 


Galaxy interacting 


Detection at 7.7 /im indicates SF 


4 


with AC 


Warped disk with SF in 3 Ho knots 


4, 11 




or tidally distorted spiral arms 


10 




SF over last 10 Myr with peak at ~ 5 Myr 


2 




Secondary peak of SF ~ 10 Myr ago 


2 




Metallicity similar to AC 


8, 12 




Optical and H I emission are lopsided 


10 



Table 3 — Continued 



Nature 


Properties 


References 




Shift between H I and stellar light in the south 


10 


HCG 31 E 


Tidal debris removed from AC 


Young SSCs (1-3 Myr and « 10 Myr) 




Material falling back to AC 


with masses 10 4 - 10 5 M© 


2 




Age of component El ~ 3 Myr 


13 




Underlying old stellar population 


2 




Two components in counter-rotation 


5 




Main component counter-rotating with respect to AC 


5 




Amplitude of rotation of main component is 25 km s _1 ) 


5 




Eastern component with recent SF 






Western component is older 


2 




Connected to AC by string of SF regions 


2 




Metallicity similar to AC and B 


12 



Table 3 — Continued 



Nature Properties References 



HCG 31 H 


Tidal debris 


Blob of size of ~ 10"in between E and F 


12 


removed from AC 


Metallicity similar to AC or higher 


12 




Higher velocity than E and F by 80 km/s 


12 




SF episode with estimated age 3-4 Myr 


12, 13 


HCG 31 F 


TDG candidate 


Several SSCs (< 4 Myr) with masses 10 4 - 10 5 M 


2 




No underlying old stellar population 


2 




Age of components Fl and F2 ~ 3 Myr 


13 




Age of component F3 ~ 6 Myr 


13 




Metallicity similar to AC 


8 




Formed in tidal debris (4 Myr) 


2 




Detected in Ha but not at 7.7/um 


4 




Two components non-connected by continuum or Ha 






within the same H I peak 


2 




No rotation, low gaseous velocity dispersion 


5 




H I kinematically detached from tidal tail 


10 



HCG 31 G 



Galaxy moving through H I mass Kinematically distinct from AC 



8 



Table 3 — Continued 



Nature Properties References 



or accreting mass from H I cloud Local maximum of common H I envelope 9, 10 

Asymmetric Ha emission 2 

SF at border facing southern tail 2 

SF over last 10 Myr with peak at ~ 5 Myr 2 

Late-type morphology with very small bulge 10 

e(a) or A+em type galaxy - with considerable ongoing star formation 13 



HCG 31 Al 

Tidal debris No rotation, low gaseous velocity dispersion 5 

in northern tidal tail of AC Detected at 7.7/un 4 



Table 3 — Continued 



Nature Properties References 



HCG 31 A2 


Tidal debris removed from AC 


Rotation pattern (amplitude 25 km s^ 1 ) 


5 


Material falling back to AC 


Counter-rotating with respect to AC 


5 




Detected at 7.7/xm 


4 


HCG 31 Q 


Galaxy 


Radial velocity compatible with group membership 


8, 10 




Local maximum of common H I envelope 


9, 10 




Spectrum dominated by A-type population 






with very strong Balmer absorption lines 


13 




Possibly post-starburst with residual star formation 


13 


HCG 31 R 


Newly identified TDG candidate 


Formed by a number of star forming regions 


13 




Radial velocity compatible with group membership 


13 




Age of component Rl ~ 2.6 Myr 


13 




Metallicity similar to AC and B 


13 




Coincides with a maximum in the H I envelope 


10, 13 




Very weak continuum 


12, 13 



Note. — The following abbreviations are used in this table: FIR (far infrared emission); SB (star- 
burst), SF (star formation), SSCs (super-star clusters), W-R (Wolf-Rayet). 

References. — (1) Johnson et al. (1999); (2) Johnson & Conti (2000); (3) Conti (1991); (4) O'Halloran 
et al. (2002); (5) Amram et al. (2004); (6) Mazzarella & Boroson (1993); (7) Yun et al. (1997); (8) Richer 
et al. (2003); (9) Williams et al. (1991); (10) Verdes-Montenegro et al. (2005); (11) Iglesias-Paramo & 
Vflchez (1997); (12) Lopez-Sanchez et al. (2004); (13) this work. 



